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ABSTRACT 

We present a comprehensive study of the abundance of carbon dioxide in exoplanetary atmospheres in hot, 
hydrogen-dominated atmospheres. We construct novel analytical models of systems in chemical equilibrium 
that include carbon monoxide, carbon dioxide, water, methane and acetylene and relate the equilibrium con¬ 
stants of the chemical reactions to temperature and pressure via the tabulated Gibbs free energies. We prove that 
such chemical systems may be described by a quintic equation for the mixing ratio of methane. By examining 
the abundances of these molecules across a broad range of temperatures (spanning equilibrium temperatures 
from 600 to 2500 K), pressures (via temperature-pressure profiles that explore albedo and opacity variations) 
and carbon-to-oxygen ratios, we conclude that carbon dioxide is subdominant compared to carbon monoxide 
and water. Atmospheric mixing does not alter this conclusion if carbon dioxide is subdominant everywhere 
in the atmosphere. Carbon dioxide and carbon monoxide may attain comparable abundances if the metallic- 
ity is greatly enhanced, but this property is negated by temperatures above 1000 K. Eor hydrogen-dominated 
atmospheres, our generic result has the implication that retrieval studies may wish to set the subdominance of 
carbon dioxide as a prior of the calculation and not let its abundance completely roam free as a fitting param¬ 
eter, because it directly affects the inferred value of the carbon-to-oxygen ratio and may produce unphysical 
conclusions. We discuss the relevance of these implications for the hot Jupiter WASP-12b and suggest that 
some of the previous results are chemically impossible. The relative abundance of carbon dioxide to acetylene 
is potentially a sensitive diagnostic of the carbon-to-oxygen ratio. 

Subject headings: planets and satellites: atmospheres - methods: analytical 


1. INTRODUCTION 

1.1. Towards robust interpretations of spectra 

There is an ongoing debate in the astrophysical com¬ 
munity on how to interpret the measured spectra of the 
atmospheres of exoplanets. The first school of thought 
uses a series of “forward models”: given a set of assump¬ 
tions about the atmosphere, one com putes forward and 


predicts a synthetic spectrum 

fe.g.. iSeager & Sasselov 

2000 

ISudarskv. Burrows & Pinto 120001 IBurrows et al. 

2007 

: iCahov. Marlev & EortnevI 

20101 lEortnev et al. 20101 

Barman et al.l 1201 It iMarlev et al. 

I2OI2I ISniegel & Burrows 

2012h. This approach has the advantage that it is grounded 


by the laws of physics and chemistry (radiative and chem¬ 
ical equilibria). It draws from a rich heritage of, and has 
enjoyed success in, interpreting brown dwarf spectra ( e.g., 

iM^lev et al.lll99d: iTsuji et ^119961: iBurrows et al.lll99% . 

The disadvantage is that Nature may outsmart our 
preconceived notions of an atmosphere. The second 
school of thought uses “atmospheric retrieval”, which 
is the attempt to shed these p reconceived notion s (e.g. , 
Madhusudhan & Seaged l2009t iBenneke & Seaged 1201^ 
Lee, Heng & IrwinI 120131 iLine et aTT l20 1 4l) . This approach 
is sound when remote sensing and in-situ measurements are 
available, such as for the Earth and the Solar System bodies, 
but its robustness is not assured when scrutinizing distant, un¬ 
resolved point sources of light. Eor e xample, workers enforce 
global energy conservation (e.g., iMadhusudhan & Seagerl 


120091) . but this does not guarantee local energy conservation 
(radiative equilibriuirQ; see Appendix). Generally, the ro¬ 
bustness of an interpretation via retrieval depends on the prior 
knowledge that one is assuming about the atmosphere—one’s 
model is only as good as the assumptions one inserts into it. 

In the present study, our goal is to study the gaseous chem¬ 
istry of hot, hydrogen-dominated atmospheres. It may not 
apply to atmospheres with rocky surfaces, where outgassing 
may produce carbon-dioxide-dominated atmospheres. By 
“hot”, we mean that water cannot condense out, at least in the 
photospheric regions observed by astronomers. Eurthermore, 
if the inert gas of the atmosphere is a mixture of hydrogen 
and helium, we assume that atmospheric escape does not act 
to remove hydrogen, leaving behind a helium-dominated at- 
mosphere that may contain more carbon diox ide than water 
(iHu & Seagedl2014lHu. Seager & Yun^l2015h . 

1.2. May we set generic priors on carbon dioxide? 

To move forward, it would be useful to uncover generic 
physical and chemical trends that one may insert into these in¬ 
version techniques as priors. In the current paper, we focus on 
the study of carbon dioxide (CO2), which is typically a minor 
carrier of carbon in exoplanetary atmospheres. Understand¬ 
ing its abundance relative to other molecules such as carbon 
monoxide (CO), water (H2O) and methane (CH4) is directly 
relevant to dec iphering the carb on-to-oxygen ratio (C/O) of 
an atmosphere (iLine et al.ll2013h . 
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C/O 


P-CH4 + nco + ncOs 
^HaO + nco + 2nco2 ' 


( 1 ) 


^ Whether radiative equilibrium is a good assumption, in three dimensions, 
is another matter altogether. But in the context of one-dimensional models, 
this is a valid discussion to have. 
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Fig. 1.— Reciprocal of normalized equilibrium coefficient associated with 
the production and destruction of carbon dioxide. 



Fig. 2.— normalized equilibrium constants as a function of temperature for 
two different values of pressure (P = 10 mbar and 1 bar). Note that has 
no pressure dependence. 

where the mixing ratios are represented by n and labeled with 
self-explanatory subscripts. Specifically, if CO is the domi¬ 
nant molecule, then we have C/O « 1. By contrast, if CO 2 
dominates, then we have C/O Ri 0.5. For example, this 
has direct consequ ences for the debate on w hether WASP- 
12b is carbon-ric h (iMadhusudhan et ^1201 Ih or carbon-poor 
jLine et al.ll2014h : IMadhusudhan et al.l (1201 ll) report CO be- 
ing do minant over CO 2 (see their Figure 2), while iLine et al.l 
( |2014|) report CO 2 being dominant over CO (see their Table 
3). 

The main conclusion of this study is that carbon dioxide 
is almost always subdominant in hot, hydrogen-dominated at¬ 
mospheres, compared to carbon monoxide and water, because 


nc02 

nco 

nc02 

%20 


^2 
nco 
K'2 ^ 




' nco- 


( 2 ) 


We will elucidate what “almost” means. At this juncture, it is 
not so important to understand what the dimensionless func- 



T(K) 



T(K) 


Fig. 3.— Benchmarking of our approximate carbon-rich (thick, dashed 
curves) versus exact carbon-poor (thin, solid curves) solutions. The approx¬ 
imate solutions assume that he, ho 1- When the carbon-to-oxygen 
ratio is low (top panel, C/O = 0.5), the solutions match perfectly; the 
acetylene abundance is computed from it not formally be¬ 

ing included in the carbon-poor solutions. When it is high (middle panel, 
C/O = 10), the abundances of methane and carbon monoxide are over- and 
under-estimated, respectively, because of the formal exclusion of acetylene; 
the estimate for acetylene, using the carbon-poor solutions, breaks down. As 
the elemental abundance of oxygen becomes large {ho = 0.1, which is a 
factor of 200 larger than the value at the solai* photosphere), our approximate 
and exact solutions begin to deviate from each other (bottom panel). 
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tion K 2 is—this will come later. It is more important to note 
that it has no dependence on pressure and that 1 jK'^ ^ 0.1-1 
for a broad range of temperatures (Figure [T]). Since we typi¬ 
cally have ri-HaO 1 and nco <C 1 in a hydrogen- or helium- 
dominated atmosphere, we may conclude that, in most situa¬ 
tions. 


»C 02 

nco ’ ^H 20 


(3) 


These trends hav e previously been elucidated b y 
MadhusudhanI (1201^ . iKop pyapu^Kasting & Zahnl^ ((20l3), 
Moses et al.l (12013*30) and IVenot et al.r(l2014li 7 in one form or 
another, using numerical calculations of chemical kinetics. 
Our present study takes a complementary approach and 
corroborates this conclusion via a set of analytical calcula¬ 
tions of equilibrium chemistry, spanning a broad range of 
temperatures, pressures, C/O values and metallicities that 
are representative of all of the currently observable and 
characterizable exoplanetary atmospheres. (Whenever it 
is relevant, we point out the limitations of our analytical 
approach.) 

Our main technical contribution is to present novel analyt¬ 
ical solutions for chemical systems with carbon dioxide, car¬ 
bon monoxide, water, methane, acetylene (C 2 H 2 ) and molec¬ 
ular hydrogen, which allow us to efficiently explore vast 
swarths of parameter space. We use tabulated values of the 
Gibbs free energy to relate the equilibrium constants to tem¬ 
perature and pressure. We then fold these calculations of at¬ 
mospheric chemistry with analytical models of temperature- 
pressure profiles to explore its effects in a wide variety of ex¬ 
oplanetary atmospheres. 


2. METHODOLOGY 

2.1. Carbon-poor (C/O < 1) chemistry 

We consider the conversion of methane to carbon 
monoxide via the net reaction (iBurrows & Shard 119991: 
iLodders & Feglevll2002llMoses et al.ll2Ql Ih. 

CH4 -f H2O ^ CO -f 3H2. (4) 


Note that the number densities of the molecules marked by 
tildes have been normalized by nB_ 2 , while those of the atoms 
will be normalized by tt-h. The former are the mixing ratios, 
while the latter are the normalized elemental abundances. 

Stoichiometric book-keeping (counting the number of 
atoms of each species) yields 


ncHi + nco + nco2 = nc, 

nii20 + nco + 2nco2 = no, (9) 

4ncH4 + 2 nH 20 + 2 nH 2 = ?^H- 


Manipulating these equations and using the equilibrium 
constants yield a cubic equation for the mixing ratio of wa¬ 
ter, 

^3ii-H20 + ^2tiH20 + CinH20 + Co = 0, (10) 

where we have defined 


C3 = — (2no - 4 nc - 1) 

A. 


C2=K' 


2 (no - nc) - 1 + ^ (n-o - 2nc) 

A, 


( 11 ) 


Cl = 2no + 4fic — 1 + 2iT' (no — n-c) 

Co = 2no- 


Unlike when carbon dioxide is excluded, we do not ob- 
tain a quadratic equation for the mixing ratio of methane 
(iHeng. Lyons & T^l2015h . We solve this cubic equation us¬ 
ing the polynomial. polyroots routine in Python. 

The rest of the mixing ratios may be obtained using 


__ 2 nc (1 -I- nH2o) _ 

ncH4 - ^ ^ ^ K'fil^o/K '2 - 4nc ’ 

nco = K'hc}iin}i 20 , 
nconH 20 


( 12 ) 


ncOa = 




This reaction alone is insufficient for modeling carbon-rich at¬ 
mospheres, because hydrocarbons such as acetylene and hy¬ 
drogen cyanide (HCN) are expected to appear at high temper¬ 
atures (lMadhusudhanll2()T^ iMoses et al.ll2013^ IVenot et H] 
l2015h . We supplement it with another net reaction for pro¬ 
ducin g carbon dioxide (ILodders & Feglevll2002t iMoses et alJ 

Honl), 

CO 2 + H 2 ^ CO -f H 2 O. (5) 

If molecular hydrogen is vastly more abundant than water, 
then the production of carbon monoxide is favored. 

The dimensional equilibrium constant of the first reaction 
is 

K' = ^ (g) 

«CH4^H20 ?T-CH4^H20 


In the low-temperature limit, we obtain fiH 20 ~ 2no and 
ncHi ~ 2fic, such that the ratio of me thane to water abun¬ 
dance s is the carbon-to-oxygen ratio (IHeng, Lyons & Tsail 
120151) . At low temperatures, water and methane each se¬ 
quester an extra atom of oxygen and carbon, respectively, at 
the expense of carbon monoxide. 

2.2. Carbon-rich (C/O >1) chemistry 

Following IHeng, Lyons & T^ (120151) , we include acety¬ 
lene as a proxy for all of the hydrocarbons that may form at 
high tem peratures and in carbon - rich situations, vi a the net 
reaction (ILodders & Feglevll2002HMoses et al,ll201 ll) , 

2 CH 4 ^C 2 H 2 -f 3 H 2 . (13) 


which may be normalized to obtain 

K' 

K' = ^. 

For the second reaction, the equilibrium constant is already 
normalized. 


(7) 


J.I _ ?T'COt^H 20 ’^C 0 ’^H 20 _ 

7^eq,2 — „ „ — - — ^ 2 - 


nn2nco2 


nco2 


( 8 ) 


Its equilibrium constant is 


^eq.3 = 


rtC2H2ttH2 _ ^C2H2^H2 


‘CH 4 


'-CH4 


which we also normalize. 


iT' =^. 


(14) 


(15) 
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TABLE 1 

Molar Gibbs free energies oe chemical species 

USED IN THIS STUDY (Pq = 1 BAR) 


T 

H 2 O 

CH 4 

CO 

CO 2 

C 2 H 2 

(K) 

(kj/mol) 

(kJ/mol) 

(kJ/mol) 

(kJ/mol) 

(kJ/mol) 

500 

-219.051 

-32.741 

-155.414 

-394.939 

197.452 

600 

-214.007 

-22.887 

-164.486 

-395.182 

191.735 

700 

-208.812 

-12.643 

-173.518 

-395.398 

186.097 

800 

-203.496 

-2.115 

-182.497 

-395.586 

180.534 

900 

-198.083 

8.616 

-191.416 

-395.748 

175.041 

1000 

-192.590 

19.492 

-200.275 

-395.886 

169.607 

1100 

-187.033 

30.472 

-209.075 

-396.001 

164.226 

1200 

-181.425 

41.524 

-217.819 

-396.098 

158.888 

1300 

-175.774 

52.626 

-226.509 

-396.177 

153.588 

1400 

-170.089 

63.761 

-235.149 

-396.240 

148.319 

1500 

-164.376 

74.918 

-243.740 

-396.288 

143.078 

1600 

-158.639 

86.088 

-252.284 

-396.323 

137.861 

1700 

-152.883 

97.265 

-260.784 

-396.344 

132.665 

1800 

-147.111 

108.445 

-269.242 

-396.353 

127.487 

1900 

-141.325 

119.624 

-277.658 

-396.349 

122.327 

2000 

-135.528 

130.802 

-286.034 

-396.333 

117.182 

2100 

-129.721 

141.975 

-294.372 

-396.304 

112.052 

2200 

-123.905 

153.144 

-302.672 

-396.262 

106.935 

2300 

-118.082 

164.308 

-310.936 

-396.209 

101.830 

2400 

-112.252 

175.467 

-319.165 

-396.142 

96.738 

2500 

-106.416 

186.622 

-327.358 

-396.062 

91.658 

2600 

-100.575 

197.771 

-335.517 

-395.969 

86.589 

2700 

-94.729 

208.916 

-343.643 

-395.862 

81.530 

2800 

-88.878 

220.058 

-351.736 

-395.742 

76.483 

2900 

-83.023 

231.196 

-359.797 

-395.609 

71.447 

3000 

-77.163 

242.332 

-367.826 

-395.461 

66.421 


Note: the molar Gibbs free energy associated with H 2 is 0 J mol ^ by 
definition. 


TABLE 2 

Molar Gibbs free energies oe the 

THREE NET REACTIONS (Pq = 1 BAR) 


T 

^Go,i 

AGo,2 

^Go,3 

(K) 

(kJ/mol) 

(kJ/mol) 

(kJ/mol) 

500 

96.378 

20.474 

262.934 

600 

72.408 

16.689 

237.509 

700 

47.937 

13.068 

211.383 

800 

23.114 

9.593 

184.764 

900 

-1.949 

6.249 

157.809 

1000 

-27.177 

3.021 

130.623 

1100 

-52.514 

-0.107 

103.282 

1200 

-77.918 

-3.146 

75.840 

1300 

-103.361 

-6.106 

48.336 

1400 

-128.821 

-8.998 

20.797 

1500 

-154.282 

-11.828 

-6.758 

1600 

-179.733 

-14.600 

-34.315 

1700 

-205.166 

-17.323 

-61.865 

1800 

-230.576 

-20.000 

-89.403 

1900 

-255.957 

-22.634 

-116.921 

2000 

-281.308 

-25.229 

-144.422 

2100 

-306.626 

-27.789 

-171.898 

2200 

-331.911 

-30.315 

-199.353 

2300 

-357.162 

-32.809 

-226.786 

2400 

-382.38 

-35.275 

-254.196 

2500 

-407.564 

-37.712 

-281.586 

2600 

-432.713 

-40.123 

-308.953 

2700 

-457.830 

-42.509 

-336.302 

2800 

-482.916 

-44.872 

-363.633 

2900 

-507.97 

-47.211 

-390.945 

3000 

-532.995 

-49.528 

-418.243 


Stoichiometric book-keeping is generalized to 

nciii + nco + PCO 2 + 2nc2H2 = ^C, 

fT-HaO + + 2nc02 = PQ, (16) 

4ncH4 + 2nH20 + 2nH2 + 2nc2H2 = n-a- 

Instead of a cubic equation for the mixing ratio of water, we 
obtain a pair of coupled quadratic equations for the mixing 
ratios of water and methane, 

^^^H 20 ^CH 4 , n ~ \ o~ ~ 

- — -h - 2nc?^H20 

+ ^icH4 + K'nii^oncUi - 4ncPCH4 - 27 tc = 0, 

2i^'nH2o”'CH4 , ^ ^2 , ~ n o- \ 

-—- 2K^nonQ^^ + n-HaO (1 - 2no) 

+ K'nn^onciii — 4noPCH4 — 2no = 0. 

The preceding pair of equations looks deceptively simple, but 
it is actually difficult to solve numerically as the mixing ra¬ 
tios span more than 20 orders of magnitude in value across 
the range of temperatures and pressures we are interested in. 
It is challenging to maintain numerical stability across such 
an enormous range of values. Taking an analytical approxi¬ 
mation is tricky, because it is difficult to judge which terms 
to drop or neglect. It turns out that we may employ an alge¬ 
braic trick if we recast the equations in (fTTl i in terms of nco. 
rather than nH 2 0 . which allows us to combine them into a 
single equation that is essentially a quadratic equation for the 
mixing ratio of CO, 

4" (4 + 2nc — 2no) — 2iT3ncH4 (1 — nc -f no) 

“ ^CH4 “ 2ncH4 (no - he) (2ncH4 + 1) = 0. 

(18) 

Notice how the mixing ratios of carbon monoxide and 
methane are no longer coupled to each other within the same 
equation—there are no “mixed” terms, unlike for each equa¬ 
tion in (fTTI l between water and methane. This property has the 
virtue that we may cleanly take the approximation he, ho 
1 and end up with a relatively simple expression for the mix¬ 
ing ratio of carbon monoxide, 

nco ~ K'nciii [ 2 ^ 6371^4 + PCH 4 + 2 (no - n-c)] • (19) 

Despite not being exact, the preceding expression retains gen¬ 
erality, because the only assumption we have made so far is 
that the elemental abundances are small compared to hydro¬ 
gen. It is less obvious how to take the he, ho -C 1 approxi¬ 
mation directly using the equations in (fTTl i. 

The expression for nco leads to an approximate expression 
for the mixing ratio of water, 

nHaO ~ 2 iT 3 ncH 4 + ficYii + 2 (no - nc), ( 20 ) 

and a quintic equation for the mixing ratio of methane, 

5 

AncH4 ~ O’ 

i=0 


( 21 ) 
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T(K) T(K) 




T(K) T(K) 


Fig. 4.— Mixing ratios of molecules as a function of temperature. We have assumed no = 5x10 ^. Top-left panel: sub-solar carbon abundance (C/0 = 0.1). 
Top-right panel: C/0 = 1. Bottom-left panel: C/0 = 10. Bottom-right panel: C/0 = 100. 


where the coefficients are 
SK'K!? 


■^5 = 


2K' 

As = —^ [1 + Sifg (no - nc)] + 2K'K'^, 

A2 = (^o - nc) + 2X3 + 

Ai = (no - nc)^ + 1 + 2 K' (no - nc ), 
.4,0 = — 2 hc- 


( 22 ) 


The mixing ratios of CO and CO 2 may be obtained using 
the second and third equations in (fT^ . respectively, while 
we have nc 2 H 2 = 4n indication that this equa¬ 

tion is correct comes from the fact that it automatically yields 
4 ch 4 ~ 2nc when all of the normalized equilibrium con¬ 
stants vanish (i.e., the low-temperature limit). Figure|3]shows 
the benchmarking of our carbon-rich versus carbon-poor so¬ 
lutions, which match perfectly when C/O < 1. 

This (approximate) quintic equation that governs the mix- 


in g ratio of methane is a nov el re sult that generalizes the work 
of iBurrows & Sha^ (119991) and iHeng. Lyons & Tsail (120151) . 
iBurrows & SharpI (I1999h excluded both acetylene and car¬ 
bon di oxide in the analytica l soluti on listed in their appendix, 
while [Heng, Lyons & Ts"^ (120151) did not consider carbon 
dioxide. 

23. Adopting the chemist’s convention for thermodynamic 
quantities 

Let the 7?.univ be the universal gas constant, p be the mean 
molecular weight, m = be the mean molecular mass, 
TOu is the atomic mass unit and fee = 7^univ/-^A be the Boltz¬ 
mann constant. If we denote the mass density by p, the vol¬ 
ume by V, the temperature by T, the number of particles by 
N and the number density by n, then the ideal gas law may 
be expressed as 

NV ■ T 

P = nknT = , (23) 

VI\a 

where Na — 6.02214129 x 10^^ mol”^ is Avogrado’s con¬ 
stant. A mole of matter contains exactly Na particles. By 
definition, TOu^a = 1 g mol~^. We write the ideal gas law 
in this way, because we wish to express the universal gas con¬ 
stant in mks units as 72.univ = 8.3144621 J mol~^ in 
order to utilize the data from the thermodynamic databases. 
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There is some confusion over the definition of the universal 
gas constant, because one may only switch between its erg 
mol”^ and erg g~^ forms, without paying a unit 

conversion penalty, when dealing with c gs units. _ 

In a departure from the formalism of iHeng. Lyons & Tsail 
(Eoil, we write G as the Gibbs free energy, rather than the 
specific Gibbs free energy. It follows that the scaling relation 
involving it becomes 

N/Na ^ N/Na ^ (^) ’ 

where Gq is the Gibbs free energy at the reference pressure 
(Po)- Typically, it is the molar Gibbs free energy at a reference 
pressure (Go = GqNa/N), and not Go, that is tabulated in 
thermodynamic databases. It has units of J mol~^. 


2.4. Equilibrium constants, Gibbs free energy, temperature 
and pressure 

We need to distinguish betwe en the different defini¬ 
tions of the equilibrium c onstant (IVisscher & Mosei 1201 It 
[Heng, Lyons & T^ |2015j) . Th e dimensionless equilibrium 
constant is (e.g., lJacobsonl2005l) 

( AGo 1 \ 

iTeq = exp - -, (25) 

\ '^univ-^ / 


where AGo.i = AGojiNA/iV and AGo,i is the change in 
Gibbs free energy going from the reactants to the products, at 
the reference pressure, associated with the first net reaction. 
iTeq should not be confused with the dimensional equilibrium 
constants we have used so far, which have been rendered di¬ 
mensionless by dividing by . 

If we focus on the first net reaction, then the dimensionless 
(Keq), dimensional (K' ) and normalized (K ') equilibrium 
constants are related by (iHeng. Lyons & T^l2015h 

K^ = nlK,^ = nl^K', (26) 

where no = Po/k-eT is the number density corresponding to 
the reference pressure and P = Thus, we obtain 


K' = 



AG, 


0,1 


TZuaivT 


(27) 


For the reaction involving carbon dioxide, we have 

^ AGo,2 


K '2 = exp - 


TZunivT 


For the reaction involving acetylene, we have 

2 




exp - 


AG, 


0,3 


TZuaivT 


(28) 


(29) 


Here, the pressure P is interpreted as the total pressure of 
the atmosphere, which is exerted mostly by H 2 in hydrogen- 
dominated atmospheres. 

Table 1 lists calculations or measurements of the molar 
Gibbs free energies of the molecules involved in the present 
study, taken from the NIST-JANAF databas^l Table 2 lists 


^ http://kinetics.nist.gov/janaf/ 


the change in molar Gibbs free energy, at the reference pres¬ 
sure of Pq = 1 bar, for all three net reactions, which allow us 
to relate K', ATg to T and P (Figure |2|i. 

2.5. Temperature-pressure profiles of model atmospheres 

To apply our calculations of equilibrium chemistry to at¬ 
mospheres, we need to know their thermal structures. To this 
en d, we employ the a nalyti cal temperature-pressure profile s 
of IHeng et d1 (120121) and Heng, Mendonca & Leei (120141), 
which generalized the pure-absorption models of iGuilloli 
(I 2 OIOI) to include non-isotropic scattering in both the visible 
and infrared range of wavelengths. These models require a 
small number of input parameters; the equilibrium tempera¬ 
ture (Peq), the interior/internal temperature (Tint), the visible 
opacity (KvIs), the infrared opacity (kir) and the Bond albedo 
(Ar). Two of these parameters are directly observable or in¬ 
ferable quantities (Tgq and Ar). In principle, the opacities 
may be the outcome of a retrieval calculation. The internal 
temperature is a quantity that is unconstrained by the obser¬ 
vations for transiting exoplanets, but we will adopt a finite 
value to include its effect. We do not explore the effects of 
collision-induced absorption, scattering by large particles in 
the infrared or Gaussian cloud decks, as these are secondary 
effects. 

It is important to note that these temperature-pressure pro¬ 
files are formal solutions of the radiative transfer equation 
and thus enforce local and global energy conservation—with 
the former being radiative equilibrium— by construction (see 
Appendix). They are not ad hoc fitting or parametric func¬ 
tions, which are commonly used in retrieval models. Rather, 
they are first-principle calculations with simplifying approx¬ 
imations taken, the most major of which is that starlight and 
thermal emission are grouped into separate wavebands. 

3. RESULTS 

3.1. Mixing ratios (molecular abundances normalized by that 
of molecular hydrogen) 

In Figure |4l we show calculations of the mixing ratios of 
carbon dioxide, carbon monoxide, water, methane and acety¬ 
lene. We set no = 5 X 10“^, which is the approximate value 
at the solar photosphere. For a sub-solar carbon abundance 
(C/O = 0.1) , water is the domin ant molecule regardless of 
temperature (iMoses et al.ll2013'bl) . Methane dominates car¬ 
bon monoxide at low te mperatures (iBurrows & Shamlll999l: 
iLodders & Feglevll2002l) and the trend reverses at high tem - 
peratures depending on the pressure (iMadhusudhanI . 
Carbon dioxide is dominated by methane at low temperatures, 
but the trend reverses in a pressure-independent manner, be¬ 
cause the number of reactants (two) versus products (two) is 
the same in the chemical reaction involving the production of 
carbon dioxide. Acetylene always remains subdominant in a 
carbon-poor environment. As we incr ease C/O to 1, 10 an d 
100, the system becomes water-poor (IMadhusudhanI . 
The trend between methane and carbon monoxide persists in 
a qualitative sense, but the transition occurs at different tem¬ 
peratures. At high temperatures, acet ylene starts to domi¬ 
nate both methane and carbon dioxid e (IMadhusudhanI l2()T^ 
iMoses et al.ll2013^ IVenot et ^l20L5l) . 

A key limitation of our approach is that graphite starts 
to condense out of the gas when the ca rbon-to-oxygen rati o 
and/or metall icity are grea t ly enha nced (IMoses et al.ll2013B) . 
For example, IMoses et akl (12013bl) haye shown that, at P = 
0.1 bar, graphite formation occurs at C/O > 0.7 and T > 700 
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Fig. 5.— Ratio of the mixing ratios of carbon dioxide to the other molecules. We have assumed no = 5 X 10 Top-left panel: sub-solar carbon abundance 
(C/0 = 0.1). Top-right panel: C/0 = 1. Bottom-left panel: C/0 = 10. Bottom-right panel: C/0 = 100. 


K when the metallicity is 300 times above solar; their Fig¬ 
ure 7 also shows that, at solar C/O, graphite forms when the 
metallicity is 1000 times solar and T < 800 K. Generally, 
graphite formation and rain ou t are expected to m aintain the 
gas-phase C/O near unity ( see iMos es et alJl 2013bl and refer¬ 
ences therein). We note that iHu & Seaged (120141) also did not 
include graphite formation in their calculations. 


3.2. Ratios of mixing ratios 

In Figure |5] we divide the mixing ratio of carbon dioxide 
by that of the other molecules. The striking thing we notice is 
that carbon dioxide is always subdominant compared to car¬ 
bon monoxide and water, across the entire range of tempera¬ 
tures considered, and for C/O = 0.1, 1, 10 and 100. Physi¬ 
cally, the relative scarcity of carbon dioxide stems from the 
fact that it both requires the production of carbon monox¬ 
ide and successfully competing with it and water for oxy¬ 
gen atoms. Furthermore, fic 02 /^CO is insensitive to pressure 
for the reason previously mentioned, unless acetylene starts to 
become dominant. By contrast, a higher pressure favors the 
product of methane over carbon dioxide as the reverse reac¬ 
tion is favored, because the number of products (four) exceeds 
the number of reactants (two)—this is a manifestation of Le 
Chatelier’s principle. 


3.3. Chemistry in a broad range of model atmospheres 

Next, we wish to compute the abundance of carbon diox¬ 
ide relative to carbon monoxide and water in model atmo¬ 
spheres and fully explore the effects of pressure variations. 
We pick C/O = 0.5 for illustration. In Figure |6] we cre¬ 
ate temperature-pressure profiles of model atmospheres with 
three different equilibrium temperatures: Teq = 600, 1200 
and 2500 K. These values are representative of a broad range 
of currently characterisable exoplanetary atmospheres from 
GJ 436b to HD 189733b to WASP-12b. For our other pa¬ 
rameters, we choose typical, plausible or illustrative values: 
g — 10^ cm s“^. Tint = 300 K, KvIs = 0.01 cm^ g“^ and 
Kir = 0.02 cm^ g~^. These choices imply that the infrared 
photosphere resides at a pressure of 

-PiR ~ — = 50 mbar. (30) 

Kir 

Starlight is most l y deposited at a pressure of 0.63p/«:vis 
(iHeng et al.ll2012t iHeng, Mendonca & Leell2014l) . For each 
value of Teq, we create three variations: pure absorption, a 
finite Bond albedo (which we set to Ar =0.5) and a temper¬ 
ature inversion (which we create by doubling the value of the 
visible opacity to K^is = 0.04 cm^ g”^)- This gives us a total 
of 9 model atmospheres. 

For each model atmosphere, we calculate hco 2 /^CO and 
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Fig. 6 .— Temperature-pressure profiles of model atmospheres with equi¬ 
librium temperatures of Teq = 600, 1200 and 2500 K, representative of 
GJ 436b, HD 189733b and WASP-12b, respectively (top panel). For each 
equilibrium temperature, we vary the Bond albedo and visible opacity of the 
atmosphere. We then calculate the abundance of carbon dioxide relative to 
carbon monoxide (middle panel) and water (bottom panel). In all panels, 
the solid curve is the fiducial model for each Teq (zero albedo, no tempera¬ 
ture inversion). The dotted curve is for Ab = 0.5, while the dashed curve 
introduces a temperature inversion by increasing the visible opacity by a fac¬ 
tor of 4. We assume a solar abundance of elements (no = 5 x 10““^ and 
C/0 = 0.5). 


nc 02 /^H 2 0- Across all 9 models (see Figure |6ll, we find that 
the mixing ratio of carbon dioxide is always subdominant, 
compared to carbon monoxide and water, by at least an or¬ 
der of magnitude—often more. Thus, we have demonstrated 
the conclusion previously stated in equation @, at least for 
elemental oxygen abundances that are solar. 


4. DISCUSSION 

4.1. Do atmospheric dynamics and mixing affect our 
conclusion? 

Our chemical-equilibrium calculations do not account 
for disequilibrium chemistry that may arise from atmo¬ 
spheric mixing. Workers constructing one-dimensional 
models of the atmosphere have traditionally used a diffu¬ 
sion coefficient, commonly termed an “eddy mixing coef- 
hcient”, to collectively mimi c advection, conve c tion, d if¬ 
fusion and turbulence fe.g., IMo ses et al l 120111 l2013allbl: 
Venot et aH 120121 l2014t ^ilndez et al.l 120141 iHu & Seageil 


2014t [Hu. Seager & Yung! T2015h . In three dimensions. 


atmospheric circulation hardly resembles diffusion with 
equator-to-pole circulation cells that penetrate down to 
^ 10 bar across the temperature range we are examin¬ 
ing dSho wmanet^ 120091: iHeng. Menou & PhillippsI 1201 1 


n l2009M Heng. Menou & Phillipps 
[Heng, Frierson & Philli ppsI [201 ipPerna. Heng & Pon 

” " " Urn' 


2012 


IKataria et al.n2013nMa\me et al.[l2014ft . In the case of GJ 
436b, the atmospheric circula tion penetrates down to only 
~ 1 bar (iLewis et al.l 120101) . (For a review of the at- 
mospheric dynamics of hot exoplanetary atmospheres, see 
Heng & ShowmanI 120151 ) Furthermore, it is often over¬ 
looked that the intense stellar heating associated with close-in, 
transiting exoplanets flattens thei r temperature-pressure pro¬ 
files and suppresses convection (iHeng, Frierson & PhillippsI 
1201 11) . Thus, the association between eddy diffusion/mixing 
and convection becomes even more tenuous. 

Nevertheless, for the sake of discussion, we invoke the 
simplest approach for including atmospheric mixing: the 
“quenching approximation”, where one locates the point in 
the atmosphere where the chemical and dynamical timescales 
equate. Above this point (i.e., at lower pressures or higher al¬ 
titudes), the chemical abundances are frozen to their quench- 
point values. A caveat is that different chemical species will 
generally have different quench points. If carbon dioxide is 
subdominant compared to carbon monoxide and water ev¬ 
erywhere in the atmosphere, then locating the quench point 
becomes irrelevant—quenching will not alter this outcome. 
Generally, the conclusion reached in equation @ is unaf¬ 
fected by atmospheric dynamics and mixing. For atmospheric 
mixing to alter our conclusion requires s ubsta ntially enhanced 
metallicities and low temperatures—in 114.31 we will quantify 
what “substantially enhanced” and “low” means. 


4.2. Does photochemistry strengthen or weaken our 
conclusion? 

Photoch emistry is expected to enhance the production of 
acetylene (iLine et al.ll201 lllMoses et al.ll201 ll) . 

2CO-F H2 ^ C 2 H 2 + 20, 

2 CH 4 ^C 2 H 2 + 3 H 2 , ^ ^ 

where the first and second net reactions proceed via the pho¬ 
tolysis of carbon monoxide (into C and O atoms, which act 
as radicals) and water (into the O atom and the OH radical), 
respectively. The first reaction is relevant only for P < 1 
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Fig. 7.— Relative molecular abundances (top-left panel) and their ratios (top-right panel) for enhanced metallicity (fio = 2.5 X 10 ^) and C/0 = 0.5, meant 
to be representative of GJ 436b. In the top-right panel, the ratio of carbon dioxide to acetylene exceeds the axis limit. Also shown are the calculations in Figure[6] 
repeated for fio = 2.5 X 10“^ and C/0 = 0.5 (bottom panels). 


/ibar, while the seco nd reaction is the d ominant one for pho- 
tospheric pressures (iMoses et al.ll201 1]) . The second net re¬ 
action is essentially the same as the one we consider in the 
present study for thermochemistry. Thus, if it is enhanced 
via photochemistry, then it can only lead to enhanced acety¬ 
lene production at the expense of carbon dioxide. Overall, the 
presence of photochemistry strengthens, rather than weakens, 
our conclusion. 

4.3. The effects of enhanced metallicity 

Our results, which are analytical except for the use of a 
numerical solver for the quintic equation, appear to be sup¬ 
ported by more sophisticated calculations in the literature 
that f ocus on somewhat narrower regimes of parameter space 
(e.g. .IMo ses et al .ll201 It i Kopparapu. Kasting & Zahnle ll2012l 
Madhusudhanl 120121: IVenot et al.l I2012L12014 FIHu & Seagerl 
2014I) . Exceptions include lLine et al.l (|201 ih . who focused on 
the hot Neptune GJ 436b and studied the effects of enhanced 
metallicities; they reported nco^/nco ^ 1 at some locations 
in their model atmo spheres with 50 x the solar metallicity. 
IMoses et al.l (12013H) demonstrated that carbon dioxide could 
be the dominant gas by mass, over a temperature range consis¬ 
tent with that found in GJ 436b, if the elemental carbon abun¬ 
dance is enhanced by four orders of magnitude over its solar 


value. lAgilndez et alJ (1201 4l) also noted that nc 02 /^H 20 in¬ 
creases dramatically with metallicity. 

We interpret an enhanced metallicity to be an increased 
value of the elemental abundance of oxygen (no) relative to 
its solar value (5 x 10“''^). In Figure |7] we repeat our cal¬ 
culations for an enhanced metallicity of 50 x solar, but we 
keep the carbon-to-oxygen ratio at its solar value (C/O = 0.5) 
as metallicity and C/O are independent quantities. We also 
focused on a lower temperature range that is representative 
of GJ 436b. We a re thus able to rep roduce the same qual¬ 
itative behavior as iLine et alJ (1201 lb : nco^/hco ^ 1 at 
T w 500-600 K. The cross-over point appears to occur 
at about 600 K. Despite the enhanced metallicity, our GJ 
436b-like calculation displays nc 02 /^H 20 ^ 1 throughout. 
We also repeated the calculations in Figure |6] and find that 
11002/^001 ^C 02 /^H 20 ^ 1 for HD 189733b- and WASP- 
12b-like atmospheres, despite the enhanced metallicity. 

Figure [8] expresses the effects of metallicity in a more gen¬ 
eral way. We examine pressures of 1 /ibar and 10 bar, be¬ 
cause these values bracket the conceivable range of pressures 
probed by the infrared observations. We examine pressures 
up to 10 bar, because this is the maximum depth to which 
atmospheric circulation penetrates for t he range of equilib- 
rium temperatures we are interested in (iPerna. Heng & Poni 
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Fig. 8 .— Mixing ratios of carbon dioxide relative to those of carbon monox¬ 
ide and water for T = 600 K (top panel) and 1000 K (bottom panel) and the 
full range of pressures that is relevant to infrared photospheric emission and 
atmospheric circulation in exoplanetary atmospheres (1 /ibar < -P < 10 
bar). We assume a solar carbon-to-oxygen ratio (C/0 = 0.5). 

I 2 OI 2 I) . We see that atmospheres that have T = 600 K always 
have nc 02 /'^H 20 < 1 and only have fico^/^CO > 1 when 
the metallicity is enhanced by more than an order of magni¬ 
tude relative to solar. By the time the temperature reaches 
T = 1000 K, we rarely have fico^/^co > 1 unless one en¬ 
tertains implausible metallicities reaching ~ 10 ^ times that of 
solar. 

A few considerations give pause to contemplating such high 
metallicities. To begin with, we note that the elemental abun¬ 
dances of carbon and oxygen in Jupiter are enhan ced by only 
about a n order of magni t ude r elative to solar (IWong et al.l 
120041) . iMiller & FortnevI (1201 Ih have shown, using interior 
structure calculations applied to mass and radius measure¬ 
ments, that gas-giant exoplanets have bulk enhancements of 
metallicities ~ 1-10 times of their parent stars. Furth ermore, 
carbon-rich stars appear to be rare (iTeske et al.ll2014h . 

4.4. Consequences and implications for retrieval studies 

Our results have important consequences for retrieval stud¬ 
ies, which is that chemical abundances cannot be allowed to 
completely roam free. For atmospheres with temperatures of 
1000 K and above, one may wish to set, as a prior in one’s re¬ 
trieval procedure, that the mixing ratio of CO 2 cannot exceed 


that of CO and H 2 O, 


nco2 

nco 


< 1 and 


nco2 

%2 0 


< 1 . 


(32) 


Retrieval outcomes that violate these criteria, for hot atmo¬ 
spheres, should be flagged as being unphysical. 

We join th e debate on the atmospher e of the hot Jupiter 
WASP-12b. iMadhusudhan et aTI (1201 Ih has previously re¬ 
ported that C/O > 1. Their Figure 1 shows four best-fit 
models with but an inspec¬ 

tion of their Figure 2 reveals that, while nco^/^co < 1 
appears to hold for their ensemble of best-fit models (which 
are colored purple), the same models als o seem to have 
nc 02 /^H 20 > 1- By contrast. iLine et al.l (l2014h reported 
that C/O < 1 for WASP-12b, but this appears to be based on 
retrievals that yield fic 02 /'^COj ^C 02 /'™H 2 0 ^ as stated 
in their Table 3 a nd dis played in their Figu re 3. Furthermore , 
IStevenson et al.l (12014l) have claimed that iLine et aH (12014h 
do not include acetylene or hydrogen cyanide in their analy¬ 
sis, which drives the retrieval calculations to “unrealistically- 
large” abundances of CO 2 . It is thus unsurprising that 
the two sets of authors reach different conclusions regard¬ 
ing the carbon-to-oxygen ratio of WASP-12b, even without 
discussing the controversy surrounding the data themselves 
(ICrossfield et al.l2012HStevenson et al.l2014l) . In other words, 
some of these reported retri evals are che mically impossible 
(but see iMoses et al.ll2013ah . To be fair, ILine et all (120141) 
did remark; “One might question how realistic solutions with 
such high abundances of CO 2 may be. Such high CO 2 
abundances are generally not thermochemically permissible 
in hig hly reducing hot-Jupiter-like atmospheres.” ILine et aH 
(120141) then performed a separate retrieyal analysis by setting 
a hard limit on the mixing ratio of CO 2 to be 10“® and found 
that it has no effect on the CH 4 and CO abundances, although 
the retrieyed H 2 O abundance is increased by 2 to 3 orders of 
m agnitude. _ 

ISteyenson et al.l (12014ft present an improyed analysis of the 
day side emission of WASP-12b, which included new data 
from the Spitzer Space Telescope. They did not enforce 
the priors in equation (l32l i. so as to allow for uncertainties 
associated with the absorption opacities and for the possi¬ 
bility of enhanced helium abundances. This permits both 
physically plausible and implausible retrieyals to be reported 
and for the reje c tion o f implausible solutions after the fact. 
ISteyenson et al.l (12014ft concluded that oxygen-rich atmo¬ 
spheres re quire implausib l y-high abu ndances of carbon diox- 
ide. Both iMadhusudhanI (120121) and ISteyenson et al.l (12014ft 
concluded that the dayside emission spectrum of WASP-12b 
fayors a carbon-ri ch interpretation. 

More recently, iKreidberg et al.l (120151) obtained a near- 
infrared transmission spectrum, which probes a series of wa¬ 
ter lines originating from the terminators of the atmosphere 
of WASP-12b, using the Wide Field Camera 3 on the Hub¬ 
ble Space Telescope. Using a thermochemically-consistent 
retrieyal analysis, they ruled out a carbon-rich scenario for the 
terminator region of WASP-12b at more than 3cr-leyel confi¬ 
dence. 

Within the context of our study, we And that the atmosphere 
of WASP-12b is too hot for enhanced metallicity to play a role 
(see Figures | 6 ] |7] and |3i in boosting the CO 2 abundance oyer 
that of CO and H 2 O. If we repeat our calculations in Figure 
0 for T = 2000-3000 K (not shown), we obtain the result 
that the metallicity needs to be enhanced by about 4 orders of 
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magnitude for WASP-12b—a statement that is independent of 
pressure, because of the high temperatures involved. Another 
consequence is that, since nc 02 /^C 2 H 2 appears to vary from 
^ 1 to ^ 1 as C/O increases, this pair of molecules may 
be u sed as a diagnost ic for the carbon-to-oxygen ratio (see 
also IVenot et al.ll2015h . although this will be challenging as 
their spectral lines will probably be dominated by those of the 
other molecules. 

Generally, past claims in the literature of nco 2 /^CO 2> 1 
and/or fico^ /^ 1 should be viewed with skepticism, 
unless a chemical mechanism (e.g., enhanced metallicity at 
low temperatures) has been elucidated. If CO 2 is reported to 
be constrained by a retrieval calculation but CO is not, then it 
should also be viewed skeptically. 


The motivation for this study originated from a discussion 
between KH and Nikku Madhusudhan during KH’s visit to 
Cambridge University in December 2014. We thank Mike 
Line, Julie Moses and Olivia Venot for constructive and help¬ 
ful comments on an earlier version of this manuscript, as well 
as the two anonymous referees for constructive and gracious 
reports. KH acknowledges financial and secretarial support 
from the Center for Space and Habitability, the PlanetS NCCR 
(National Center of Competence in Research) framework, the 
Universities of Bern and Zurich and the Swiss-based MERAC 
Foundation. KH is grateful to Claudia for encouragement and 
support. 
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APPENDIX 


NOTES ON GLOBAL VERSUS LOCAL ENERGY CONSERVATION IN ONE-DIMENSIONAL ATMOSPHERES 


In the limit where advection, “PdV work” and thermal conduction may be ne glected, the first law of thermodynamics, which 
is a statement of energy conservation, states that (Heng, Mendonca & Leell2()14h 


^ ^ 

dt cp dm 


(Al) 


where T is the temperature, t denotes the time, cp is the specific heat capacity at constant pressure and m is the column mass. 
The quantity is the net flux, which is the difference in flux entering and exiting each layer of the model atmosphere, integrated 
over all wavelengths. In a hydrostatic atmosphere, the column mass and pressure (P) are straightforwardly related by P = rhg, 
where g is the surface gravity. 

Since opacities are temperature-dependent and temperatures depend on the opacities, one needs to iterate between the m in a 
model atmosphere. With each iteration, one may compute the fluxes in each layer and thus the net flux. By using equation (lAll i. 
one may then update the temperature-pressure profile during each iteration. The model atmosphere reaches a steady state, where 
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the temperature in each layer ceases to change, when the gradient of the net flux is zero. This is known as radiative equilibrium. 
It is a statement of local energy conservation, since it applies to each layer of the atmosphere. 

Global energy conservation involves considering the energy budget of the atmosphere and ensuring that the total amount of 
energy entering and exiting it must be conserved. By taking the first moment of the r adiative transfer equation, one may show 
that the energy content, per unit mass and time and integrated over all wavelengths, is dHeng, Mendonca & LeellMTl) 

dT 

Q = ^=ti3J- 4KpasBr4, (A2) 

am 

where kj is the absorption mean opacity, Kp is the P lanc k mean opacity, J is the total intensity integrated over all wavelengths 
and (TSB is the Stefan-Boltzmann constant. Equation (IA2b still applies in a layer-wise fashion, with the first term (after the second 
equality) accounting for scattered light and also flux entering and exiting the layer from adjacent layers. The second term accounts 
for the thermal emission of e ach layer. Within the context of the method of moments, we emphasize that no approximations have 
been taken on equation (IA2b . 

To account for global energy conservation, we need to integrate the model atmosphere over column mass. 


/ Q dm = Q (m, oo), 
J m 

where we have defined, for any arbitrary parr of column masses rhi and m 2 , 

/•m2 

Q{rhi,m 2 ) = Q dm. 


(A3) 


(A4) 


We now have the mathematical machinery to describe local versus global energy conservation. Local energy conservation or 
radiative equilibrium occurs when Q = 0. This in turn implies that <5(0, 00 ) = 0, as shown by equation ( IA3b . But enforcing 
global energy conservation alone ((3(0, 00 ) = 0) does not guarantee radiative equilibrium (Q = 0). Mathematically, the gradient 
of the wavelength-integrated net flux may conspire to be non-zero at multiple locations, only to cancel out when summed over 
the entire atmosphere, such that <3(0, 00 ) = 0 but Q ^ 0. 

In other words, radiative equilibrium is a necessary and sufficient condition for global energy conservation. However, global 
energy conservation is an insufficient condition for local energy conservation. 





